Enhancement of Macroscopic Quantum Tunneling by Landau-Zener Transitions 
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Motivated by recent realizations of qubits with a readout by macroscopic quantum tunneling in a 
Josephson junction, we study the problem of barrier penetration in presence of coupling to a spin-^ 
system. It is shown that when the diabatic potentials for fixed spin intersect in the barrier region, 
Landau-Zener transitions lead to an enhancement of the tunneling rate. The effect of these spin 
flips in imaginary time is in qualitative agreement with experimental observations. 

PACS numbers; 73.40.Gk,03.67.Lx,85.25.Cp 



X3 

o 

> 

m 

O 
m 
o 



o 



X 



Macroscopic Quantum Tunneling (MQT) in Joseph- 
son systems Q has been studied in detail both experi- 
mentally and theoretically JJ in the eighties but has 
gained renewed interest very recently 0, U since the ex- 
ponential dependence of the MQT rate on parameters 
allows for a high fidelity readout of qubits based on su- 
perconducting circuits. In this context a new variety of 
quantum tunneling problem arises, namely, barrier pene- 
tration in presence of coupling to a spin-^ describing the 
qubit. For fixed spin the phase of the Josephson junction 
tunnels through a diabatic potential barrier, and the ex- 
ponentially large difference between the MQT rates cor- 
responding to the barrier potentials for the two spin di- 
rections is exploited in the readout to determine the spin 
state. However, there is an interesting parameter range 
where the diabatic potentials intersect in the barrier re- 
gion so that Landau-Zener (LZ) transitions may arise. 
In contrast to the well-studied problem of LZ transitions 
in real time [ij, here, within the standard semiclassical 
approach to quantum tunneling, these transitions occur 
in "imaginary" time. It is shown that they cause an en- 
hancement of the tunneling rate in the appropriate pa- 
rameter range in agreement with experimental observa- 
tions 0- 

To motivate the Hamiltonian underlying this work, we 
briefly consider the quantronium circuit consisting of a 
small superconducting electrode, the island, carrying an 
excess number of N Cooper pairs. This island is attached 
to a superconducting loop via two Josephson junctions 
with Josephson energy Ej/2, and it can be biased by a 
voltage source U through a gate capacitance Cg inducing 
the dimensionless gate charge Ng = CgU /2e. For suffi- 
ciently large Coulomb charging energy Ec = (2e)^/2C, 
where C is the island capacitance, and for Ng near i 
only two charge states TV = 0, 1 are relevant, and the is- 
land corresponds to an effective spin-i system that can 
serve as a qubit manipulated via the gate. For read- 
out purposes the superconducting loop is interrupted by 
a Josephson junction with Josephson energy E'j that is 
shunted by a large capacitance C" and can be biased by a 



current lb- This circuit is described by the Hamiltonian 



n = Ec{N - Ngf - Ej cos 



cos((5) 
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- - E'j co8{e) ~ -^9 



with the conjugate observables [N, 6] = —i for the qubit 
and [Q,9] = —2ie for the readout junction, respectively. 
(j) denotes the external magnetic flux through the super- 
conducting loop in units of the flux quantum h/2e. 

Within the subspace spanned by the eigenvectors 
|0),|1) of the operator N, the qubit can be described 
by Pauli matrices ai. Further, measuring all energies in 
units of Ej, we obtain the dimensionless Hamiltonian 



H = ecTz — j cos 



2m 



cos{0)~ih9 (1) 



with the dimensionless parameters e = {Ec / Ej){Ng — ^), 
j = Ej/Ej, and it = hIb/2eEj. The variables pe = 
Q/2e and 9, with the commutator [pe,^] = —i, may be 
viewed as dimensionless momentum and coordinate of a 
particle with dimensionless mass ni — C'E'j/Ae^. We do 
not discuss here the manipulations of the qubit done for 
ib = 0, but address the readout when the bias current 
is increased to a value slightly below the dimensionless 
critical current if, = 1. The problem at hand then is 
tunneling of this ficticious particle through a barrier of 
the potential cos(0) — ib9 in presence of the interaction 
with the spin. 

In view of the large mass m the coordinate 9 is almost 
a classical variable. When the kinetic energy pg/2m is 
neglected, the Hamiltonian ^ can easily be diagonalized 
with the eigenvalues 



A±(«fc, 



- cos(6i) - ib9 ± y/e^+f cos2[(6' + ?!))/2] 



that determine two adiabatic potential surfaces. At zero 
temperature and for vanishing bias current the system 
will approach the minimum of the lower surface X-{ib = 
0, 9) which for small e and j lies close to 9 — 0. As usual 
in MQT experiments, the switching of the bias current 
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form to a value close to ic = 1 is slow compared to the 
characteristic time scales of the circuit. When the system 
follows the bias current adiabatically, the particle lies at 
finite bias if, near the minimum 9-(ib) of the adiabatic 
potential X-{ib,9). This state serves as the initial state 
for the calculation of the tunneling rate. 

For this initial state it is natural to use the spin eigen- 
vectors |6'_, +), — ) associated with the eigenvalues 
X±{ib,6-) as a basis for a matrix representation of the 
Hamiltonian We then find 



H = 



v 




(2) 



where 



v± {0) = - cos{9) -ibe± (^^f^+v^ + Ko [v{e) - Vo]^ 

(3) 

are now two diabatic surfaces corresponding to the two 
spin orientations +), — ), and 



A{0) = EKo 



Vn 



is the 9 dependent coupling between them. For conve- 
nience, we introduced 

F(0)=jcosp + <^)/2], 

and Vq = y{9-), as well as 



Apparently, for 6 = 6- the Hamiltonian ||2J) is diagonal. 
It further becomes diagonal in the limits e ^ or j ^ 0. 

Depending on the external flux (/) and e/j the two di- 
abatic potentials Q may intersect. As can be seen from 
Fig. 1 this is always the case for near ^ and e small 
compared to j. Near such a crossing point 9^,, which is 
determined by 



1/2 



the diabatic potential surfaces are strongly coupled by 
the off-diagonal element of the Hamiltonian ||2J). To quan- 
tify the strength of this coupling we may introduce the 
parameter 



9(0) 



V+i9)^V^i9) 



(4) 



which diverges at the crossing point 6**. A LZ region 
where spin-flips may occur can then be deflned by the 




FIG. 1: Diabatic potential surfaces Q in dependence on e/j 
and (j) for a bias current of it = 0.96. In the region denoted by 
I the potentials do not intersect as shown in the upper inset 
displaying the potentials for e/j = 0.4, = 0.5. In region II 
the potentials cross in the barrier region as shown in the lower 
inset for e/j = 0.1, (j) ~ 1-5. The phases 0- in the metastable 
initial state before tunneling and 0, at the crossing point are 
also indicated. The dotted line refers to Fig. 3. 



condition g{9) > 1. We assume that this region is re- 
stricted to the vicinity of 9^. In fact, when the bias cur- 
rent is sufficiently far from the critical current to allow 
for several states in the metastable minimum, the LZ re- 
gion turns out to be narrow except near the boundary 
between the regions I and II in Fig. 1. This crossover 
region will be discussed further below. 

To determine the^tunneling rate we employ the 
"bounce technique" 8] which relates the rate F essen- 
tially to an imaginary time trajectory in the inverted po- 
tential, the so-called bounce. This method is equivalent 
to WKB and starts out form the partition function of the 
metastable system 



Z = Tr{|6'_,-)(6i_ 



-t3H 



} 



which has to be evaluated in the semiclassical limit for 
f3 —> oo. Within the path integral representation this 
takes the form 



Z 



V[9] e-^ 

'13 



V / dS2n--- / dsiA[0(s2„)]---A[0(si)] 
n=l Jo Jo 



X exp 



E / dT{vAe]-v+[9]) 



where the path sum runs over all orbits with period (3 
switching 2n times between V_ and V+ at times si < 
S2 < . . . < S2n- The Euclidian action on is 



S-= dT 
Jo 
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In the semiclassical hmit Z decomposes into ~ 
Ztjj~\-Zq~\-Zi. Here Z^ is the partition function of the weU 
which is obtained by summing over paths in the vicinity 
of the trivial trajectory 6(t) = 0_ sitting at the well min- 
imum. Zq is the contribution of paths in the vicinity of 
the standard non-flip bounce trajectory in In region 
II we also have to take into account the contribution Z2 
from paths that flip when the bounce traverses the LZ 
region. Trajectories with four and more spin-flips can 
be neglected away form the boundary between regions I 
and II. Both the bounce and spin-flip bounce are sad- 
dlepoint trajectories with an unstable fluctuation mode 
which after an analytical continuation yield imaginary 
and, compared to Z^^, exponentially small contributions 
that determine the rate. Following standard procedures, 
we obtain for the dimensionless rate, in units of E'j /Ti^ 



F = hm -3(Zo + Z2)/Z^ 



Fn 



(5) 



To evaluate this explicitly, we first note that within 
the barrier region the potentials V± can very accurately 
be approximated by cubic polynomials. It is convenient 
to introduce for each diabatic potential the frequency at 
the well bottom 



and a scaled distance 



v^{e±)/m , 



0±)/iol-0±) 



between the well bottom and the "exit point" 0^. where 
V±iS±) = V±{0±). Further, V,^ is the barrier height with 
respect to the minimum V±{0±). 

Since we are interested in the limit /3 — > 00, it is natural 
to look for solutions in the time interval s £ [—(3/2, (3/2]. 
The simple bounce trajectory in the inverted potential 
-VI (6*) then reads (see Fig. 2) 

^bouncc(s) ^0- + {0l - 0_)/cOsh2(w_ s/2). 

This trajectory dominates the non-flip contribution Zq 
that has been evaluated previously yielding the well- 
known MQT rate in the absence of damping 



Fo = 6V6c^-n"/7r exp(-^^ 

UJ— 



(6) 



which determines the rate in region I of Fig. 1 . 

In region II the semiclassical trajectory may switch to 
the potential surface V+{0) in the LZ region. The lowest 
order flip contribution reads 

Z2 = [ V[0] e'^-^"^ r ds2 r dsiA[0{s2)]A[0{si)] 
J Jo Jo 



X exp 



dT{V.[0] - V+[0]) 



To determine Z2 we proceed as follows: First, the ac- 
tion for the flip bounce is calculated for arbitrary flipping 
times si < S2. Due to energy conservation and the pe- 
riodic boundary condition one finds that the restriction 
0{si) — 0{s2) applies. For a path running in the inter- 
val [—(3/2, (3/2] this means that the flips have to occur 
symmetrically around s — 0. Second, as a function of si 
the action has a minimum at an optimal flipping time s| 
determined by 



Pe^{Vm-V-[0]} 







(7) 



where t^[0] is the time the bounce spends on the V+ 
surface. As one might have guessed, Eq. {Tj) yields 
V+[0] = V-[0], so that the optimal flips occur at the 
intersection point 0* of the diabatic potentials. Then 
0* = 0{sl) and S2 - si = := <+. [Other solu- 

tions of Eq. (0 with pe{sl) — mean that flips occur 
either in the well or at the turning point. In both cases 
the orbit has no energy to run on V+ and one regains 
the simple bounce action on Further, solutions with 
dt^[0]/d0 — corresponds to a maximum of the action.] 
This way the trajectory at the saddlepoint of the action 
is obtained as 



0fi 



.p(^) = h{[.s[ - f) 



^mp(^) + My 



Here the step functions h{-) select the time segments 
spend on the two potential surfaces, where 



t4 



0^4- A 



F(ip[m). 



(8) 



The parameter A = [P'{xo)]^^'^ is determined by the slope 
P' of the polynomial 



P{x) = — — — p, with p 



v+{0. 



V-(0-] 



(9) 



at its zero where P{xo) — 0, and F((p|m) is the ellip- 
tic integral of the first kind with modulus m — 1/2 + 
P"{xo)/8X^ and angle cos{ip) = [A^ - {xq - xt)]/[X^ + 
(xq - xt)]. Further, 0^ip{s) = 6'bounco[s - sign(s)a] de- 
scribes the segments of the flip bounce on the surface V- 
where it coincides with the simple bounce apart from a 
phase 



a 



■ Arccosh( V-T* ) 



Finally, the segment running on V+ follows as 

,cn{XsLL}j^[fh) — 1 



.+p(s) ^0+ + {0l - 0+) 



xa + A^ 



cn(Asa;+|?T!,) + 1 



where cn(-) is a Jacobi function (see Fig. 2) 
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FIG. 2: Non-flip bounce (thin line) and flip bounce (thick line) 
vs. time. The inset shows the inverted potential surfaces —V- 
(thin line) and — V+ (thick line), see text for details. 



For the action of this flip bounce one finds 



18V^ 

5 w_ 

5uj+ 



^J{l-x-)3 (^3y^(l-a;,-)-5j +2 
2xo - +9p 



2A 



F((^|?fi) + 2AE(</?|m) 



2 * \ X'' + xq - x2 



P"{< 



Here, E(-) is the elliptic integral of the second kind. 

Next, the action is expanded around S<, with respect 
to variations of the flipping times up to second or- 
der. For this purpose it is convenient to introduce sum, 
u = si + §2, and difference, v = S2 — si, times, respec- 
tively. The Gaussian factor describing deviations from 
s*,S2 then depends only on v, i.e. only on the relative 
position of the flips but not on the absolute position of 
the bounce in time. Hence, like the conventional non-flip 
bounce, the flip bounce has one zero mode (corresponding 
to the integration over u), while the integration over v is 
weighted with a Gaussian factor exp[— 17^/2 (u — 1+[6'*])^]. 
The frequency 



m 



dt+{e) 



de 



is proportional to the second derivative of the action at 
s*, Here, t^{d) for arbitrary 9 is given by the expres- 
sion ((SJ where the parameters \(9) and fh{0) are gained 
by replacing in ® p p[e) = p + [V-{0) - V+{9)]/V^+ 
and the angle (p{9) follows from cos((p) = [A^ — {xq — 
x+)]/[X^ + (xo - x+)] with x+i9) = {9- 9+)/{9% - 9+). 

Finally, fluctuations in 9 are calculated in the usual 
way j3| from a ratio of two determinants where the zero 
eigenvalue is omitted and replaced by the proper zero 
mode normalization factor while the unstable mode is 
accounted for by an analytic continuation leading to the 
imaginary part. 

Eventually, one obtains in the limit (3 oo for the 
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FIG. 3; Rate enhancement vs. magnetic flux along the dotted 
line in Fig. 1. Arrows indicate the boundary between regions 
I and II. 



contribution r2 to the rate ||SJ) 



Fa =6i 



A (61, 



rfc 



V2 



This can be combined with ® to yield the central result 
of this paper, namely, the decay rate in region II 



F = rn 



with S- 



1 



A(6i, 



— erfc 
2 



nit+{9,y 


e(s-s.)| 


[ V2 J 
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Vf^ /lo^. The corresponding rate enhance- 
ment is shown in Fig. 3. Apparently, there is a pro- 
nounced exponential increase of the total rate due to 
spin flips in the LZ range along the bounce. This is in 
agreement with experimental flndings of the peak current 
variation with magnetic flux in the quantronium device 

As one approaches the boundary between I and II the 
LZ region grows and in a narrow boundary layer the pa- 
rameter g{9) in is larger than 1 in the entire barrier 
range. Then, multi-spin flips can occur anywhere along 
the bounce and also during the switching on of the bias 
current. However, since the exponential factors of Fq and 
F2 coincide at the boundary, this breakdown of the nona- 
diabatic approach utilized here, essentially reduces to a 
prefactor effect smoothing the transition between the re- 
sults in region I and II. A more interesting extension of 
the present work would consider initial states where due 
to manipulations via the charge gate the system moves 
on the upper adiabatic potential surface A-|_(ib,0). This 
will be addressed in future work. 
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